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A formula for the two-loop infrared singularities of dimensionally regularized QCD scattering am- 
plitudes with an arbitrary number of massive and massless legs is derived. The singularities are 
obtained from the solution of a renormalization-group equation, and factorization constraints on 
the relevant anomalous-dimension matrix are analyzed. The simplicity of the structure of the ma- 
trix relevant for massless partons does not carry over to the case with massive legs, where starting 
at two-loop order new color and momentum structures arise, which are not of the color-dipole form. 
The resulting two-loop three-parton correlations can be expressed in terms of two functions, for 
which some general properties are derived. This explains observations recently made by Mitov et 
al. in terms of symmetry arguments. 



I. INTRODUCTION 

In the past few years, remarkable progress has been 
achieved in the understanding of the infrared (IR) singu- 
larities of massless scattering amplitudes in non-abelian 
gauge theories, which are characterized by an intricate 
interplay of soft and collinear dynamics. While factor- 
ization proofs guarantee the absence of IR divergences in 
inclusive observables [lj, an all-order formula for the IR 
singularities of QCD amplitudes has been lacking for a 
long time. An important step toward this goal was made 
by Catani [2J , who correctly predicted the singularities of 
two- loop scattering amplitudes apart from the 1/e pole 
term. An interesting alternative approach to the problem 
of IR singularities was developed in [3[, where the au- 
thors exploited the factorization properties of scattering 
amplitudes along with IR evolution equations to recover 
Catani's result at two-loop order and relate the coeffi- 
cient of the 1/e pole term to a soft anomalous-dimension 
matrix. This matrix was later calculated at two-loop or- 
der in 0, Q . Its color structure was found to be of the 
same form as at one-loop order and it was shown that 
the more complicated color structure of the 1/e-term in 
Catani's formalism Q was an artifact of his subtraction 
scheme. 

In recent work Q , we have shown that the IR singular- 
ities of on-shell amplitudes in massless QCD are in one- 
to-one correspondence to the ultra-violet (UV) poles of 
operator matrix elements in soft-collinear effective theory 
(SCET) US EH- They can be subtracted by means of a 
multiplicative renormalization factor, whose structure is 
constrained by the renormalization group. We have ar- 
gued that the simplicity of the corresponding anomalous- 
dimension matrix holds not only at one- and two-loop 
order, but is in fact an exact result of perturbation the- 
ory. A first test of this prediction at three-loop order was 
performed in [111 ]. Detailed theoretical arguments sup- 
porting our conjecture were presented in where we 
used constraints derived from soft-collinear factorization, 
the non-abelian exponentiation theorem [13l . |14| , and the 
behavior of scattering amplitudes in two-parton collinear 
limits [HI to show that the anomalous-dimension matrix 



retains its simple form at least to three-loop order, with 
the possible exception of a single color structure multi- 
plying a function of conformal ratios depending on the 
momenta of four external partons, which vanishes in all 
collinear limits. Some of these arguments were devel- 
oped independently in We also showed that higher 
Casimir invariants, which in particular would lead to a 
violation of Casimir scaling of the cusp anomalous dimen- 
sion [3, E3j EH El ) do not appear at four- loop order. 

It is interesting and relevant for many physical appli- 
cations to consider generalizations of these results valid 
in the case of massive partons. The IR singularities of 
one-loop scattering amplitudes containing massive par- 
tons were obtained some time ago in [20(. We will re- 
produce their formula as a special case of our analysis. 
In the limit where the parton masses are small compared 
with the typical momentum transfer among the partons, 
they serve as regulators for collinear singularities. In 
[2ll . I22I [23j factorization theorems for this limit were pro- 
posed, which allow one to derive massive amplitudes from 
massless ones. This technique has been used to derive the 
massive e + e~ — > e + e~ scattering amplitude [H, [24[ and 
the virtual corrections to heavy-quark production in the 
limit where all kinematic invariants are larger than the 
heavy-quark masses [25], [26[ . 

First steps toward solving the problem of finding the 
IR divergences of generic two-loop scattering processes 
with massive partons, without restricting oneself to the 
limit of small masses, have recently been taken by Mi- 
tov et al. [27]]. Interestingly, these authors find that the 
simplicity of the anomalous-dimension matrix governing 
the IR poles does not persist at two-loop order in the 
massive case. Specifically, they point out that there ex- 
ist non-vanishing singularities in Feynman graphs con- 
necting three different partons. They evaluate the cor- 
responding contributions numerically and point out that 
they vanish for some special kinematic configurations. 

In this paper we extend our analysis in [l2| to the 
general case of gauge-theory amplitudes with arbitrary 
numbers of massive and massless partons. We show that 
the structure of the terms found in [13] follows from sim- 
ple symmetry arguments, such as soft-collinear factoriza- 
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tion in SCET and non-abelian exponentiation. The fact 
that in the presence of massive partons the low-energy 
effective theory knows about the 4-velocities of heavy 
particles in addition to the light-like directions of mass- 
less partons gives rise to additional color and kinemati- 
cal structures, which are absent in the case of massless 
partons. As a result, the structure of IR poles becomes 
increasingly more complicated in higher orders of pertur- 
bation theory. We present, for the first time, a general 
formula for the IR singularities of dimensionally regular- 
ized two-loop scattering amplitudes with arbitrary num- 
bers of massive and massless partons and arbitrary values 
of the parton masses. It generalizes the one-loop result 
of [20]. For amplitudes with n > 4 partons, our result 
contains two new functions with certain symmetry prop- 
erties, for which at present no analytical expressions are 
available. 

We begin by considering the case where the parton 
masses are of the same magnitude as the typical mo- 
mentum transfer between the partons. In this case the 
appropriate low-energy effective theory is a combination 
of SCET and heavy-quark effective theory (HQET) [H], 
which is applicable since the relative velocities of the 
heavy partons are of 0(1) in this case. With the gen- 
eral result at hand, we then explore the limit where the 
parton masses are taken to be much smaller than the 
hard momentum transfers between the partons. We con- 
clude that if the two new functions did not vanish in 
this limit, then the QCD factorization formula of [221.1231] 
would need to be modified. 



II. SOFT-COLLINEAR FACTORIZATION 

We denote by \M n (e, {p}, {m})}, with {p} = 
{pi,...,p n } and {m} = {mi,...,m n }, a UV- 
renormalized, on-shell n-parton scattering amplitude 
with IR singularities regularized in d = 4 — 2e dimen- 
sions. This quantity is a function of the Lorentz invari- 
ants Sij = 2o ij pi ■ pj + iO and p\ = mf , where the sign 



factor <Ji 



-1 if the momenta pi and pj are both in- 



coming or outgoing, and — — 1 otherwise. We assume 
that all of these invariants are of the same order and refer 
to them as hard scales. For massive partons (rrii ^ 0), 
we define 4-velocities Vi — pi/mi, whose components are 
of 0(1). We have vf = 1 and define the abbreviations 
Wij = —a a Vj ■ Vj — iO. We use the color-space formal- 
ism of [29|, [3(| , in which n-particle amplitudes are treated 
as n-dimensional vectors in color space. Ti is the color 
generator associated with the i-th parton and acts as a 
matrix on its color index. The product T t ■ Tj = T°- 
is summed over a. Generators associated with different 
particles trivially commute, Ti ■ Tj = Tj ■ Ti for i ^ j, 
while Tf = Ci is given in terms of the quadratic Casimir 
operator of the corresponding color representation, i.e., 
C q = Cq = Cf for quarks and C g = Ca for gluons. 

We have shown in 0, [12] that the IR poles of such am- 
plitudes can be removed by a multiplicative renormaliza- 



tion factor Z~ l (e, {p}, {m}, /i), which acts as a matrix 
on the color indices of the partons. This quantity obeys 
the renormalization-group equation 

Z- 1 ^ Z(e, M, {m}, M ) = -T(M, {m}, /i) , (1) 

where T is the anomalous-dimension matrix of effective- 
theory operators built out of collinear SCET fields for the 
massless partons and soft HQET fields for the massive 
ones. The formal solution of this equation is 



Z{e,{p},{m},p) = Pexp 



r^r({ £ },{m},//) 

Jut 



m 

where the path-ordering symbol P means that matrices 
are ordered from left to right according to decreasing val- 
ues of n' . The ^-factor appearing in the renormalization 
of effective-theory operators describes the IR behavior of 
on-shell amplitudes, because these amplitudes are closely 
related to the bare Wilson coefficients of the correspond- 
ing operators. This connection is discussed in detail in 
[12l ]. Compared to the massless case studied there, we en- 
counter one complication: Since virtual corrections due 
to heavy quarks are integrated out in the effective theory, 
the strong coupling constant entering the ^-factor in the 
low-energy theory is defined in a theory with massless 
quark flavors only, while the massive amplitudes in QCD 
also receive contributions from heavy-quark loops. The 
Z-factor we obtain from the effective theory describes 
the IR singularities of massive QCD amplitudes after the 
coupling constant is matched onto the effective theory 
with massless flavors. The corresponding matching rela- 
tion will be given below. 

The interactions between collinear and soft fields can 
be decoupled by means of a field redefinition Q, after 
which soft interactions manifest themselves as interac- 
tions between a set of light-like and time-like soft Wilson 
lines representing the massless and massive particles, re- 
spectively. Generalizing the discussion of [I4], the rele- 
vant soft operator in the present case is 



S({n},{E},A0 



■S Vn \0), (3) 



where partons 1, . . . , k are massless, and the remaining 
n — k partons are massive. 

From now on, we label the massive partons by capi- 
tal indices /, J, . . . and the massless ones by lower-case 
indices i, j, . . . . The anomalous-dimension matrix of the 
effective-theory operators can be written as a sum over 
soft and collinear contributions [Hj], 

r(M, {™}, p) = r s ({^}, fi) + J2 ito, m) , (4) 

i 

where T s is the soft anomalous-dimension matrix govern- 
ing the UV poles of the Wilson- line operator in ([3]). The 
collinear contributions T* only arise for massless partons 
and are diagonal in color space. Note that color conser- 
vation implies the relation 
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Tr = 



(5) 
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when acting on color-singlet states such as color- 
conserving scattering amplitudes. In intermediate steps 
in the calculation of the anomalous-dimension matrix one 
needs to regularize IR divergences in the effective theory, 
for instance by taking the massless partons slightly off 
their mass shell, (— pf) > 0. Following [L2\ . we intro- 
duce the notation Li — ln[/i 2 /(— pf )] for the associated 
collinear logarithms, which need to cancel in the final re- 
sult ([4]). The soft anomalous-dimension matrix T s is, in 
the most general case, a function of the cusp angles flij, 
Pij, and flu formed by the Wilson lines belonging to 
different pairs of massless or massive partons. With the 
IR regulator as specified above, the relations expressing 
these cusp angles in terms of the hard momentum trans- 
fers and particle masses read 1 



cases where at least one massless parton is involved has 
been explained in [12, EH • 

It is instructive to see how the dependence on the IR 
regulators disappears in (|U), when we combine the ex- 
pressions in ([7]) and ([8]) and express the cusp angles in 
terms of hard momentum transfers and masses as well as 
collinear logarithms. We note that (after the cusp angles 
have been eliminated) 



dT s 



2 — parton 



dLi 



■Ti 



Cj 7cusp( a s) — q l 



/cusp 



K) 



(9) 
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Pi j = In ■ 



(-Pj) 



L j - In 



flu — arccosh(wjj) = arccoshf 



rriip 

- s i 3 
-sij 



(6) 



Hence, the sum of all contributions is indeed independent 
of the IR regulators. Note that this requirement fixes the 
relative strength of the terms proportional to Ti ■ Tj and 
Tj ■ Tj in ©. From (0} we then obtain 



V 2m/m 



r(W,{m},M)! 

_ ^2 Ti ' Tj 



2 — parton 



7cus P (a s ) In ■ 



The anomalous dimensions appearing on the right-hand 
side of (@| are functions of the cusp angles flij and the 
collinear logarithms L^, while that on the left-hand side 
only depends on the hard scales and m,. 



III. TWO-PARTON CORRELATIONS 

We begin by considering the one- and two-particle 
terms in the anomalous-dimension matrix (U) . They can 
be written as [3l| 



E Tl 2 rj 7cus p^ /j ' a ^ + 

Ti ■ T 3 7cusp(a s ) In 

1,3 



E 

i 

E 



(10) 



-Sij 



where j l = 7* + 7*. 

The anomalous-dimension coefficients 7cusp(o! s ) and 
"f l (a s ) (for i = q,g) have been determined to three- 
loop order in [12] by considering the case of the mass- 
less quark and gluon form factors. For example, the 
one- and two-loop coefficients in the perturbative series 



P c (Li,/i) = -Ci^ cusp (a s ) Li +7*(a s ), (7) 7cusp (a s ) = J2 n 7™ sp (fj)" +1 are 



and 



Pa ~ E Tl ? T ' 7cu Sp (/3/j, a s ) 



r «({^}.M)| 2 _ parton 

E Tj ■ Tj 
2 7cus] 

~ E T i ' T J 7cuspK) pjj + E 7l(a s ) + E 7 I (« S ) , 
1,3 <■ 1 

(8) 

where the notation (ii, refers to unordered tuples of 
distinct parton indices. The various coefficients are func- 
tions of the renormalized coupling a s = a s (p) and, in 
the case of 7 CU sp(/3, a s ), of a cusp angle fl. The fact that 
only a linear dependence on the cusp angles is allowed in 



cusp 

7o 



7i 
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(11) 

In QCD only quarks can be massive, and the first two 
coefficients in the expansion of 7^ can be extracted 
by matching our result with the known form of the 
anomalous dimension of heavy-light currents J n i in SCET 
" SI. We obtain 



Tjhib^jM) = -Cf 7cu Sp (a s ) In 



2v • p 



•YK), (12) 



where we assume that the heavy-quark with velocity v 
is incoming and the light quark with momentum p is 
outgoing. The sum 7' = r y q + 7^ was first obtained at 
two- loop order in [33[ . Using this result leads to the one- 
and two-loop coefficients 



1 Strictly speaking, in the effective theory only the large light- 
cone components of the collinear momenta appear in the scalar 
products pi ■ pj and vj ■ pj , see ll jf . 



7o Q 



-2C F , 



-'■' = 0-r,('^-f -4&)+y C F T F n f . 



(13) 
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Since this anomalous dimension is connected with the 
soft Wilson-line operator in ([3]), its color structures are 
constrained by the non-abelian exponentiation theorem 
[l3l Il4| . This explains the absence of a C\ term, and it 
implies that up to 0{o? s ) the corresponding anomalous 
dimension for a massive color-octet particle is given by 
Ca/Cf times the coefficients shown above. 

It remains to determine the velocity-dependent func- 
tion 7cus P (/3, ce s ) in (fT0|) . This can be accomplished by 
applying our general result to the special case of just two 
heavy quarks, where it should reproduce the velocity- 



dependent anomalous dimension of heavy-quark currents 
in HQET [H, GUI- It follows that 

r Jhh (v,v', a s ) =C Flcusp (P,a s )+2 1 Q(a s ) (14) 

with cosh/3 = v ■ v', where v is the incoming and v 1 
the outgoing quark's velocity. This anomalous dimension 
was calculated at two- loop order in [35l . [36j , and it has 
recently been recomputed and expressed in compact form 



371 ]. Using the result Q13p for 7^, we obtain from the 



latter paper 



7cus P C3,Q!s) = 7cus P (o! s )/3coth/3+ ^ -j coth 2 (3 

+ coth /3 



Li 3 (e-^)+pU 2 (e-^)-C 3 + ^(3+{ 

6 3 
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Li 2 (e- 2/5 ) - 2/3 ln(l - e~ 2p ) - — (1 + (3) - /3 2 - '— 



3 



(15) 



For small cusp angle the result can be expanded in even 
powers of j3, 



7cusp(/?, u s ) 



/3 2 

7cusp(d! s ) ( 1 + — 



(16) 



C A [- 
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Ca-l 
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Note that the leading terms in the expansion of 
7cusp(/3> Q-s) around small cusp angle are equal to — 2jq, 
so that the full anomalous dimension r,/ hh (v,v',a s ) van- 
ishes in the limit v ■ v' = 1. In the limit of large cusp 
angle one finds 



7cus P (/3, a„) = 7cus P (a s ) P + ■ ■ ■ 



(17) 



where the dots represent terms that vanish for /3 — > 00. 
Note that no constant terms remain in this limit. 

The .Z-factor associated with a renormalization-group 
equation such as ([T]), in which the anomalous dimension 
is linear in ln/x, was derived in Q. To two-loop order it 
reads 



4tt V4e 2 2e 



(18) 



3/3oT' , r;-4/3or , r! 



where 

r'(« s ) 



d 



16e 3 



T({p},{m},fx) 



16e 2 4e 



Tcusp 



(19) 

and we have expanded T = T fj + Ti (fj) 2 + . . . . Expo- 
nentiating the result (|18p yields the two-loop expression 



for Z, which encodes the IR singularities of the massive 
QCD scattering amplitudes. More precisely, we have 



Z- x (a s )\M n (e,{p},{m}))\ 



<<x s 



finite (20) 



for e — > 0. The quantity a s denotes the strong coupling 
constant in the effective theory, which is obtained after 
integrating out the heavy quark flavors. It is obtained 
from the coupling constant ap CD of full QCD via the 
decoupling relation a^ CD = £a: s . To first order in a s , the 
matching factor appropriate for heavy-quark flavors 
reads 1381 



rih 
i—1 



r(e) 



(21) 



We have checked that at one-loop order our results (flT 
and ([18]) reproduce the IR pole terms obtained in [2fJ. 
We have also confirmed that they correctly describe the 
IR singularities of the two-loop massive quark form factor 
calculated in [H, 0, HO] . 

Our result pTJ]) has been derived under the assumption 
of hard parton masses, in which case the appropriate low- 
energy effective theory is well known. However, being 
an exact result in perturbative QCD, this formula can 
also be applied in cases where some or all of the parton 
masses are much smaller than the momentum transfers 
between the partons, m/mj <C s/j|. This limit is rel- 
evant to many processes, such as Bhabha scattering or 
heavy-quark production at the LHC. It is the limit in 
which all cusp angles, also those involving massive par- 
tons, become large, and ftu = ln(— su/mimj) up to 
power-suppressed terms. In this case our general result 
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fLOf implies that 

r(M,{m^o}, M )| 



E 



2 — parton 



r(W,{o} )M ) 



Cj 7cu S p(as) In h 7 / ("s) - l l {a s ) 

mi 



(22) 



where T({p}, {0}, /i) is the anomalous-dimension matrix 
in the massless case, whose conjectured all-order form 
is given by the terms shown in the second line in fTO]) 
7, 12]. In the equation above, 7* is the massless single- 
parton anomalous dimension belonging to parton /. In 
QCD only quarks can be massive, and this result can be 
rewritten as a sum over heavy-quark anomalous dimen- 
sions 



FQ(mQ,(x) = C F j cusp (a s ) In — +7 y (a s ) -j q (a s ), 

m Q 

(23) 

where the one- and two-loop coefficients of the constant 
terms are 



7o 
7? 



it 



(24) 
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The factor Zq associated with (f23|) . which is obtained 
after substituting the anomalous dimension Tq into 
the general relation (fT5|) , is compatible with the re- 
sults of [H, |2j|. Specifically, we find that the product 
z -2 z {m\o} is fimt6; where the quantity ^{™|o} was de _ 

fined in [22j as the ratio of the massive to the massless 
quark form factor in the limit where the quark mass tends 
to zero, and without including heavy fermion loops. Note 
that our derivation assumed that the massive partons are 
heavy enough to be integrated out in the low-energy the- 
ory using (|2~Tjl . If this is not the case, then the treatment 
of the heavy-flavor contribution is more complicated [23| . 



IV. THREE-PARTON CORRELATIONS 

It was observed in [27j that in the case with massive 
partons the anomalous-dimension matrix ^ has a more 
complicated structure than in the massless case, and that 
at two-loop order non-abelian diagrams connecting three 
partons give rise to non-vanishing contributions. The ad- 
ditional terms were found to vanish if two of the three 
partons are massless, 2 or if any pair of the three kine- 
matic invariants formed out of the parton momenta are 



2 It is noted in |27j that this observation has been made indepen- 
dently by Einan Gardi. 




FIG. 1: Graphical representation of the two three-particle 
terms in the anomalous-dimension matrix (|27p . Double lines 
represent massive partons, single lines show massless ones. 



equal. We will now show that these observations have a 
simple explanation. 

Adapting the diagrammatic analysis of our paper [l2| 
to the case with non-zero parton masses, we find that 
additional structures arise from two-loop order on, the 
reason being that the 4- velocities of the massive partons 
are known to both the full and the effective theories. In 
HQET the velocities appear as labels on the effective 
heavy-quark fields [H, E2] • In the full theory, they are 
simply given by Vi = pi/rrii. While for massless partons 
the rewriting from hard to soft variables always intro- 
duces collinear logarithms, this is not true for massive 
partons, as shown in ([6|). At two- loop order, the non- 
abelian exponentiation theorem then allows additional 
structures involving three partons. They are absent in 
the massless case, because it is impossible to form a to- 
tally anti-symmetric function of three cusp angles fyj, 
Pjki Pki that is independent of collinear logarithms upon 
the substitution shown in the first line in ([6]) [12J. This 
would violate soft-collinear factorization. However, with 
massive partons this argument no longer applies. In 
fact, in principle the soft anomalous-dimension matrix 
can contain the structures 



(I,J,K) 

+ i f abc E E T i t j T * F 2 ifiu - fa ,Pik) (25) 
(i, j) k 

+ *r bc E E T i T ! T k ^ {fin , fa, m . 

The function F± must be totally anti-symmetric in its 
arguments, while F2 (F3) must be anti-symmetric in the 
last (first) two arguments. Soft-collinear factorization 
enforces that after elimination of the cusp angles using 
the result (j2"5|) must be independent of collinear log- 
arithms. This in turn requires that 



FaifiijjPikjPjk) 



f2(0IJ,0Jk 

0, 



(26) 



6 



where fzix,])) must be an odd function of y. Note that 
for F 3 to be independent of collinear logarithms it should 
be a function of the combination + — 0jk)-> but 
this is symmetric in j, k and so vanishes when contracted 
with the anti-symmetric color structure. Hence only the 
two possibilities illustrated in Figure[T]remain, and we are 
led to the following additional structures in the complete 
anomalous-dimension matrix in (|U): 

r ({£M^}>M)| 3 _ parton 

= lf abc J- TfT^F^j^jK^Kl) (27) 
(I,J,K) 

-vjk vj ■ Pk - 



*/ a " C EE T ? T J T kh (frj, In- 

(I, J) k 



-Oik VI ■ Pk 



This is the most general form possible at two-loop order, 
and hence the sum of (|10[> and (|2~7| gives the complete 
answer for the anomalous-dimension matrix at 0{a 2 s ). 
We note that the color factor in (f27|) is non-zero only if 
there are at least four partons involved in the scattering 
process. For three partons, charge conservation, i.e. the 
fact that Ti + T 2 + T 3 = 0, implies that 



abc 



rj-ia rj-ib rj-ic 



Since the color matrices entering the two-particle and 
three-particle terms do not commute, the path ordering 
in becomes important in the massive case. However, 
since the three-parton color structures first enter at two- 
loop order, this complication arises only at 0(a^.). 

It is interesting to observe that the two structures in 
([27)) are consistent with the constraints following from 
the behavior of scattering amplitudes in the limit where 
two (or more) massless partons become collinear. To see 
this, consider the case where the first two light partons 
become collinear, such that Pj = zP and p2 = (1 — z)P 
with P 2 — -> 0. As shown in pjj, the anomalous dimen- 
sion of the matrix of splitting amplitudes Sp({pi,p2j-, £i) 
is given by the difference of F({pi,p 2 ,P3, ■ ■ ■ }, A 1 ) 
and r({P,p3, . . . }, {to}, ^)|t p =Ti+t 2 , and it must be in- 
dependent of the momenta and colors of the remaining 
partons. Indeed, the three-parton term (|2T|) does not 
contribute to this difference, because it is at most lin- 
ear in the color generators of massless partons and in- 
variant under rescalings of their momenta. It would be 
interesting to explore whether additional constraints on 
the anomalous dimension arise from considering quasi- 
collinear limits [2(| EH], a generalization of the usual 
collinear limit to the case where some of the partons in- 
volved in the splitting are massive. 

Our result p7|) explains the observation made in [27| . 
that the three-parton correlations do not arise when two 
or three of the involved partons are massless (i.e., the fact 
that -F3 = 0). Moreover, when any pair of kinematic in- 
variants are equal, these terms vanish as well. For F\ this 
follows from the fact that it is an anti-symmetric function 
in all of its arguments. For f% it follows since for vi = Vj 
both arguments of the function vanish, but f 2 must be 



odd in its second argument. We can thus reproduce and 
understand all of the observations made in [27J based on 
symmetry properties and without an explicit two-loop 
calculation. This also demonstrates that these observa- 
tions will continue to hold in higher orders of perturba- 
tion theory. On the other hand, due to reasons similar 
to the ones outlined above, it is clear that when massive 
particles are present more and more complicated color 
and momentum structures will arise in higher orders of 
perturbation theory. The arguments of [12( show that all 
of these structures are of non-abelian origin and involve 
three or more partons. 

An integral representation for the function F\ can, in 
principle, be extracted from expressions derived in [271 ] . 
These authors have shown by numerical evaluation that 
Fi / for generic values of its arguments. It is an 
open question, however, how the functions Fi and f 2 
behave in the limit of large cusp angles, corresponding 
to the case where m/raj <C I s / J I- If either one of the 
two functions does not vanish in this limit, then the fac- 
torization theorem for massive amplitudes proposed in 
[22I [HI would have to be modified to account for the 
non-factorizable three-particle terms derived from (|27[) . 
For example, in heavy-quark production processes such 
as qq — > ti or gg —* ti only the second term in [)27p 
contributes. If the asymptotic behavior of the coefficient 
function were f2(x,y) ~ xy, then in the limit s 3> m 2 
the three-parton term would contribute a 1/e IR pole at 
two- loop order proportional to ln(t/tt) ln(— s/m 2 ), which 
would be incompatible with a simple factorization for- 
mula. The agreement of the numerical results presented 
in [3] for the qq — > ti scattering amplitude with the 
predictions obtained in [25| using the factorization theo- 
rem for massive partons provides some evidence that the 
three-parton contributions do indeed vanish in the limit 
of large cusp angles, but one should check this with an 
explicit calculation of F\ and f 2 . 



V. SUMMARY AND OUTLOOK 

We have derived the general form of infrared singu- 
larities of two-loop QCD amplitudes involving arbitrary 
numbers of massive and massless partons. These singu- 
larities can be absorbed into a multiplicative .Z- factor, 
which fulfills a renormalization-group equation. In the 
purely massless case, the associated anomalous dimen- 
sion involves only color-dipole correlations. This was first 
observed by explicit calculation at the two-loop level in 
0, Q. In a recent paper, we have argued that factor- 
ization constraints suggest that this property persists to 
all orders and have analyzed the three-loop case in detail 
[T2I ]. More recently, Mitov et al. [27| have shown that 
the anomalous-dimension matrix for massive partons in- 
volves color correlations among three partons, in contrast 
to the massless case. They have evaluated these contri- 
butions numerically in Euclidean space and pointed out 
that they vanish for certain kinematical configurations. 
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In the present paper, we have derived the structure 
of the two-loop anomalous-dimension matrix relevant for 
the general case involving both massive and massless par- 
tons. The factorization constraints can be studied using 
a combination of soft-collinear effective theory for the 
massless partons and heavy-quark effective theory for the 
massive ones. We have first derived the form of the two- 
parton correlations and extracted the necessary anoma- 
lous dimensions from known results for heavy-to-heavy 
and heavy-to-light currents in the effective theory. The 
factorization constraints are weaker for terms involving 
massive legs, since both the full and the effective theories 
know about the four- velocities of the massive partons. In 
particular, the constraints do not exclude three-parton 
correlations if at least two of the partons involved are 
massive. We have shown that two such structures ap- 
pear and that their symmetry properties imply that they 
vanish when two four-velocities of the involved massive 
partons become equal, which explains the findings of [27| . 
Starting from their results, one should derive explicit rep- 
resentations for these two functions in Minkowski space, 
preferably in analytic form. 

We have briefly discussed the limit in which the parti- 
cle masses are small compared to the momentum trans- 
fers. In this case, the two-parton contribution factors into 
a hard part, depending on the large momentum transfers, 
and a sum over collinear contributions depending on the 
parton masses. It would be interesting and important to 
study this case in more detail. In particular, one should 
check whether the three-parton correlations vanish in this 
limit, since their presence would violate the factorization 
theorem proposed in [22|, [23| . 
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Note added 

Soon after this paper was published, the two-loop ex- 
pressions for the functions F± and ]i describing the three- 
parton correlations in (|2"T)) have been calculated in closed 
analytic form [45|, |4||. In light of these results, some 
statements made in the last two paragraphs of Section ITVl 
need to be revised: 

1. The conclusion drawn in [27j and in the present 
work, that the three-parton terms vanish whenever the 
four- velocities of two massive partons coincide, is false. 
Due to the presence of Coulomb singularities the limit 
Vi — y vj is singular. The three-parton terms do not van- 
ish but diverge logarithmically in this limit. 

2. A formal argument presented in [27| . suggesting that 
the three-parton terms only receive contributions from 
the three-gluon graphs shown in Figure [1] is invalid. 
Additional contributions with the same color structure 
arise from two-loop planar and one-loop counterterm di- 
agrams. The calculations presented in that paper are 
therefore not sufficient to extract an integral representa- 
tion for the function F±. 

We also note that the explicit expressions for the func- 
tions Fi and /2 vanish like (mjmj / su) 2 in the limit 
of small parton masses. The factorization theorem for 
amplitudes with massive partons proposed in [22l . |23| is 
therefore not invalidated by their presence. 
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